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Abstract 

m ■ 

Using a nonperturbative method based on asymptotic behaviour 
of Wilson loops we calculate masses of glueballs and corresponding 
\ Regge-trajectories. The method contains no fitting parameters and 

^ ; the mass scale is fixed by the meson Regge slope. Theoretical pre- 

r~| ■ dictions for lowest glueball states are in excellent agreement with lat- 

Oh, tice results. The leading glueball trajectory and its relation to the 

Q_i| Pomeron is discussed in details. Important role of mixing between 

glueball and qq trajectories is emphasized. 

.^J ; Glueballs are among the most intriguing objects both in experiment [1] 

^ ■ and lattice data [2,3]. While experimental situation is not yet settled, lat- 

tice simulations yield an overall consistent picture of lowest (< AGeV) mass 
spectrum. The mass scale and level ordering of the resulting glueball spec- 
trastrongly differ from those of meson spectra, yielding a unique information 
about the nonperturbative (NP) structure of the gluonic vacuum. Another 
important feature of glueball spectrum is that the leading glueball Regge 
trajectory is closely connected with the Pomeron singularity, which deter- 
mines asymptotic behaviour of high-energy diffractive processes - the object 
of active studies during last decade. 

In this paper we will study the problem of spectra of glueballs and of 
the Pomeron singularity using the method of Wilson-loop path integrals de- 
veloped in papers [4,5,6]. The method is based on the assumption of an 
area law for Wilson loops at large distances in QCD, which is equivalent 
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to the condition of confinement of quarks and gluons. It is shown below 
that long distance dynamics is defined by the string tension only and yields 
spin-averaged glueball masses, while spin splittings are due to mostly pertur- 
bative exchanges and relatively small. Our predictions for masses of lowest 
spin-averaged glueball states in units of yfa are in an excellent agreement 
with results of recent lattice calculations [2,3] e.g. the calculated mass ratios 
M(G)/M(0~ + ) are within 5% on average from the lattice data of [2]. Spin- 
orbit and spin-spin splittings are calculated and found in a good agreement 
with lattice data also (within la deviation on average). 

To find the Pomeron singularity, we calculate the leading glueball Regge- 
trajectory in the positive t region and extrapolate it to the scattering re- 
gion of t < 0. The importance of mixing among this trajectory and qq- 
trajectories (/, /') is emphasized. A calculation of these mixing effects leads 
to the Pomeron trajectory with op(0) > 1 (op(0) = 1.1 4- 1.2) in accord 
with experimental observations [7]. An interesting interplay of 3 vacuum 
trajectories in the region t > is observed. 

A previous study of glueball spectrum in a similar method was done in 
[8]. In this letter we improve results of [8] in several ways. First, the analysis 
of perturbative gluon exchanges (PGE) shows that PGE do not sum up to 
the adjoint Coulomb interaction, but rather to the BFKL ladder, where loop 
correction strongly reduces the sum, so that PGE can be disregarded in the 
first approximation. Secondly, we use the rotating string Hamiltonian of 
[5,6] which yields correct string slope l/2na (as compared to l/8<7 in [8]). 
Finally the mixing of glueball and /, /' trajectories is taken into account, 
which drastically changes the resulting Pomeron. 

Following ref . [4] , we separate gluonic fields into nonperturbative back- 
ground _B M and perturbative gluons a M , = B^ + a^ and consider two-gluon 
glueballs, described by the Green's functions 

G^Ax,y\x',y>) = (T m G^(x,x')G vul (y,y')T out ) + perm (1) 

where r^ n ' out ^ are operators projecting given quantum numbers and is 
the gluon Green's function of field a M in the background field B^, namely 

G^(x,y) = (x\(-D\ u - 2igF^)- 1 \y) (2) 

where = <9 M — igB^, is the field strength of the field B^ in the ad- 
joint representation, and averaging over background B^ is implied by angular 
brackets. 
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Expression (1) can be written as the path integral (see refs. [4,8]) 



G^ vl {x,y\x',y') = const J ds J ds'DzDz'e- K ~ K (W F ) (3) 
where K = | J s (^) 2 <ir, K' is the same with primed z, r, s, and 

(HV) = trP B P F (exp{tg f B^du^ + 2ig f Fdr + 2i 5 f Fdr 1 }) (4) 

jo Jo 

Here P B , Pp are ordering operators of the color matrices and F respec- 
tively. The terms with F generate spin-dependent nonperturbative contri- 
butions, which are calculable and relatively small, and we shall treat those 
terms perturbatively. 

Neglecting F's as a first approximation, one arrives at the Wilson loop 
in the adjoint representation, for which one can use the minimal area law, 
confirmed by numerous lattice data [9] up to the distance of the order of 
lfm, 

(W adj ) = Z exp(-a adj S min ) (5) 

where we have included in Z self-energy and nonasymptotic corrections, since 
(5) is valid for large loops with size r, T ^> T g , where T g is the gluon corre- 
lation length. 

Applying now the general method of [5] to the Green's function (3), one 
defines the Hamiltonian H for 2 gluons without spin corrections. 

ffo=4+**)+ L{L+1) 



l*(t) ^ 1 r*\pL{t) + 2fi(0-\) 2 v(P,t)d0\ 

+ f'^T^-y + l- [\(f3,t)dP (6) 



/o 2v(f3,t) 2 Jo 

Here L is the angular momentum in the cm. system of two gluons, L = 
[r x p], r is the distance between gluons, and p,p r are momentum and its 
radial projection respectively; t and f3 are the cm. time and distance along 
the string - two coordinates on the world sheet of the string. 

Positive auxiliary functions ji{t) and v(/3,t) are to be found from the 
extremum condition [5]. Their extremal values are equal to the effective 
gluon energy (fi) and energy density of the adjoint string (is). The physical 
meaning of those is better illustrated by limiting cases of the free gluon, 
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where (//) = (y/jP) and of the free string without quarks, where (v) 



^(l-^-i) 2 )- 1 / 2 . 
For the case L = the extremization over \x and v yields a simple answer, 
coinciding with the Hamiltonian of the relativistic potential model 

H = 2\Jtf + a adj r (7) 

With the replacement of the operators fj,(t),v(/3,t) (which by extremization 
are expressed through operators p, r) by c -numbers, to be found from ex- 
tremization of eigenvalues of H , one obtains a simple form, used in [8], 

P 2 

H = h no + (T adj r (8) 

The eigenvalues of (8) are about 5% higher than those of H . 

The value of <j a dj in (8) can be found from the string tension Of und of qq 
system, multiplying it by |, as it follows from Casimir scaling observed on 
the lattice [9]. Taking experimental Regge slope for mesons a' = 0.89 GeV 2 
one obtains af und = 0.18 GeV 2 and a adj 0.40 GeV 2 . 

In what follows the parameter \i and its optimal value /i , which is com- 
puted from eigenvalues of (8) and proportional to ^JCadj, play very important 
role. The way they enter spin corrections (see below) and magnetic moments 
shows that Ho is an effective (constituent) gluon mass (or constituent quark 
mass in the equation for the qq system) . 

Note that our lowest "constituent gluon mass" fj, (n = L = 0) = 0.528 GeV 
(for o f = 0.18 GeV) is not far from the values introduced in the potential 
models, the drastic difference is that our /i is calculable and depends on n, L 
and grows for higher states. 

The mass spectrum of H , Eq. (6), and was studied in [5,11]. With a 
good accuracy it can be described by a very simple formula 

M 2 

= L + 2n r + ci (9) 



where L is the orbital momentum, n r -radial quantum number and c\ is 
a constant ~ 1.55. It describes an infinite set of linear Regge-trajectories 
shifted by 2n r from the leading one (n r = 0). The only difference between 
light quarks and gluons is the value of a, which determines the mass scale. 
A similar spectrum was obtained independently by numerical quantization 
of qq sysytem in [12]. 
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The lowest glueball state with L = 0, n r = has M = 2.01 GeV. It 
corresponds to a degenerate ++ and 2 ++ state. 

In order to compare our results with the corresponding lattice calculations 
[2,3] it is convenient to consider the quantity M/^/oJ, which is not sensitive 
to the choice of string tension a Q. We also introduce the spin averaged mass 
M, which for L = 0, n r = states is defined as M = §(M(0++) + 2M(2++)), 
and in a similar way for higher states. This definition takes into account 
the structure of spin-splitting terms (see below in (10)), so that M can be 
compared with the eigenvalues Mo of spinless Hamiltonian (6). 

The comparison of our predictions for spin averaged masses of the lowest 
glueball states with corresponding lattice results is given in Table 1. For the 
average mass with L = 2,n r = lattice results are limited to the state 3 ++ . 

It follows from Table 1 that the spin-averaged masses obtained from 
purely confining force with relativistic kinematics for valence gluons are in a 
good correspondence with lattice data, which implies that PGE effects are 
suppressed. We come back to this point in what follows. 

Now we shall consider spin splittings of glueball masses and shall treat 
spin effects as a small perturbation. Such a treatment is justified aposteriori 
by our results and by lattice data, which demonstrate that spin splittings for 
glueball states (apart from 2 ++ — ++ ) amount to less than 10-15% of the 
total mass. 

To proceed one should choose among two possibilities of spin state de- 
scription of bound gluons. One can insist on transversality condition ( valid 
for a free gluon) also for gluons in the glueball, as it is done in the potential 
model of [10]. Instead in our formalism the bound gluon becomes massive 
due to the string ( with the effective energy /io not equal to the gluon momen- 
tum |ki|) and has 3 spin polarizations, similarly to W^^Z® bosons, which 
get their mass from Higgs condensate. For more discussion see [11]. 

Then two-gluon mass operator can be written as 

M = M (n, L) + SLM SL + S {1) S (2) M S5 + M T , (10) 

where S = S*- 1 -* + S^ 2 \ Mo is the eigenvalue of the Hamiltonian H = Hq + 
AH, pert , and H is given in (7) (or its approximation in (8)), while AH pert is 
due to PGE. 

1 Note that the values 07 ~ 0.21 — 0.23 GeV 2 used in lattice calculations differs by 
about 20% from the "experimental" value at = 0.18 GeV 2 . 
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To obtain three other terms in (10) one should consider averaging of the 
operators F in the exponent of (4) and take into account that 



2^ = 2(S«B« + S (1) E( 1 )) 



(11) 



and similarly for the term in the integral / Fdr', with the replacement of 
indices 1 — > 2. Here gluon spin operators are introduced, e.g. 



{^rn)ik — i&miki h & — 1) 2, 3, (S^)j4 — i&im (12) 



Now the spin splittings are obtained from (4) where the proper time 
(r, t') is replaced by the real (Euclidean) time (t, t') via the relation [4-6] 
dr = T^y, where //(£) is the same operator as in the Hamiltonian (6) with 
its extremal value fio playing the role of the gluon constituent mass ( and 
explicitly calculated from (6) or (8)). As a result the spin-dependent part 
of the Hamiltonian(lO) can be written in the form similar to that of Eichten 
and Feinberg [12] (see [11] for more details) 



where S = ' + S' , AV contains higher cumulant contributions which can 
be estimated of the order of 10% of the main term in (4) [4]. 

The functions Vi(r) are the same as for heavy quarkonia [13] and the 
NP part of spin splittings can be in this way expressed through the field 
correlators D(x), D 1 (x) [14], measured on the lattice [15]. The resulting 
estimates show that the only appreciable NP spin-splitting is due to Vi, Vz 
(Thomas precession) and was computed numerically as in [13]. 

Perturbative contributions to spin splittings can be calculated in the same 
way using results of refs [10]. Combining all corrections and values of M 
from Table 1 we obtain glueball masses compared with lattice data in Table 
2 for 07 = 0.23 GeV 2 . 

The general feature of spin-dependent contribution AH S is that it dies 
out with the growing orbital or radial number. This feature is well supported 
by the lattice data in Table 2. 




+ ^(3(S (1) n)(S (2) n) - S (1) S (2) )\/ 3 (r) + AV 



(13) 
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One can see in Table 2 that magnitude of our spin-spin and spin-orbit 
splittings ( proportional to (5^(r)) and (1/r 3 ) respectively) is in good agree- 
ment with lattice data (our splittings are on average some 30% higher than 
lattice spin splittings defined with rather large uncertainties); inclusion of 
PGE in the form of adjoint Coulomb would increase splittings several times 
[8] and destroy the agreement, yielding another phenomenological argument 
against adjoint Coulomb form of PGE. 

In many analytic calculations of glueball masses it is postulated that there 
is a Coulomb-type interaction between valence gluons, which differs from the 
qq case by the same | factor as in eq. (9). If this would be true it would 
lead to a drastic decrease of glueball masses calculated above due to Coulomb 
attraction, especially for L=0 (for a s = 0.3 this mass drops down by 0.5 GeV 
[8]). This clearly contradicts lattice calculations and casts a doubt on the 
validity of the assumption about presence of Coulomb interaction (without 
higher loop corrections). 

A careful study of general derivation of Coulomb-like small distance con- 
tribution in this formalism shows [11] that contrary to the quark case for 
valence gluons perturbative gluon exchanges do not sum up to the Coulomb 
potential, but rather to the BFKL ladder [16]. 

In order to estimate effects of small distance contributions we shall use the 
analysis of these effects on gluonic Regge-trajectories not from the glueballs 
mass spectra at positive t, but for t = 0. Extensive calculations of the 
gluonic Pomeron trajectory intercept have been carried out in the leading 
log approximation (LLA) [16] and a s cocorrections were calculated recently 
[17]. 

An intercept of the leading Regge pole is strongly modified by a s correc- 
tions 

A = a P (0) - 1 = a s — /n2(l - Ca s ) 

The coefficient C is large (« 6.5) and a s correction substantially decreases 
A compared to LLA result. 

The coefficient C is rather large (^ 6.5) and the a s correction strongly 
reduces A. Its value depends on the renormalization scheme and scale for 
a s . In the "physical" (BLM) scheme values of A are in the region 0.15 -j- 0.17 
[18]. We can estimate mass-shift of the lowest glueball state due to PGE 
effects using this result and assuming that the slope a' P ~ 0.4 GeV 2 of 
eq.(ll) is not be strongly modified by perturbative effects. Thus one can 
expect that characteristic shift due to perturbative effects in M 2 [L = 0, n r = 
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0),5M 2 « A/a' adj « (0.38 -=- 0.48) Gel/ 2 . This corresponds to the shift 
in M(L = 0,n r = 0),<JM « 5M 2 /2M « 0.1 Gel/. This shift is much 
smaller than naive estimate using adjoint Coulomb interaction and gives 
some justification to the neglect of PGE in our calculations. It should be 
noted that this is only a rough estimate of the perturbative effects. 

The three-gluon system may be considered in the same way, as it was 
done for the two-gluon glueballs. The 3g Green's function is obtained 
as the background-averaged product of 3 one-gluon Green's function, in full 
analogy with (1). Assuming large N c limit for simplicity and neglecting spin 
splittings one arrives at the path integral (cf equation (3)) 

3 poo 

G (3fl) = const 17/ dsiDz®e- Ki - aSi (14) 
t=i Jo 

where a = o-f und , since every gluon is connected by a fundamental string 
with each of its neighbours. 

Using as before method of refs. [4-6] one obtains the following Hamilto- 
nian (we assume symmetric solution with equal /^(t) = = 1,2,3 (no 

orbital excitations was assumed as in (8)) 

= + f + *(\ri- r 2 | + |r 2 - r 3 |) (15) 

Here p , are Jacobi momenta for the 3-body system of (q, q, g). 

To simplify treatment further, we shall consider fjj clS ct constant to be 
found from the extremum of eigenvalues, as in (8), which in that case provided 
some 5% increase in eigenvalues and we expect the same situation in this case. 

The eigenvalues of were found using the hyperspherical method in- 
troduced in [19]. A reliable estimate of M gives for the lowest state (and 
a f = 0.18 GeV 2 ) the value M = 3.22 GeV (3.62 for a = 0.23 GeV 2 ). Per- 
turbative hyperfine interaction increases the mass of the 3 spin splittings 
compared to this value by 0.29 GeV. The large value for the lowest three- 
gluon state is in an agreement with lattice calculations (see Table 2). For 
more details concerning 3g states see [11]. 

Let us consider now the problem of the Pomeron Regge trajectory in 
more details, taking into account both nonperturbative and perturbative 
contributions to the Pomeron dynamics. 

The large distance, nonperturbative contribution gives according to Eq.(12) 
for the leading glueball trajectory (n r = ) 

a P {t) = -ci + a' P t + 2 (16) 
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In eq.(16) we took into account spins of "constituent" gluons, but ne- 
glected a small nonperturbative spin-spin interactions. For an intercept of 
this trajectory we obtain ap(0) ~ 0.5, which is substantially below the value 
found from analysis of high-energy interactions [7] a p (0) — 1.1 1.2. 

The perturbative (BFKL) contribution leads to a shift (increase) of the 
Pomeron intercept by ~ 0.2 as it was explained above. The resulting a~p (0) ~ 
0.7 is still far from experiment. 

The most important nonperturbative source, which can lead to an in- 
crease of the Pomeron intercept is, in our opinion, the quark-gluon mixing or 
account of quark- loops in the gluon "medium". In the 1/N -expansion the 
effect is proportional to Nf/N c , where Nf is the number of light flavours 

In the leading approximation of the l/iV c -expansion there are 3 Regge- 
trajectories with vacuum quantum numbers, - gg-planar trajectories (a/ 
made of uu and d d quarks, af made of ss-quarks) and pure gluonic trajectory 
- ac- The transitions between quarks and gluons ~ will lead to a mixing 
of these trajectories. Note that for a realistic case of G, f and /'-trajectories 
(Fig.l) all 3-trajectories before mixing are close to each other in the small 
t region. Trajectory of gluonium crosses planar / and /'-trajectories in the 
positive t region (t < 1 GeV 2 ). In this region mixing between trajectories is 
essential even for small coupling matrix gik(t). Lacking calculation of these 
effects in QCD we will consider them in a semi-phenomenological manner. 

Denoting by a.i the bare /, /' and G-trajectories and introducing the 
mixing matrix gik(t)(i,k = 1,2,3) we obtain the following equation for de- 
termination of resulting trajectories after mixing [11] 

f ~ f Yl S < + jQ2 - 9 2 ik) - ai«2«3 + Yl ~ 2 9l29l3923 = 

For realistic values of gik{t) (for details see [11]) the resulting trajectories are 
shown in Fig. 1 by solid lines. The pomeron intercept is shifted to the values 
«p(0) > 1. For t > lGeV 2 the Pomeron trajectory is very close to the planar 
/-trajectory, while the second and third vacuum trajectories - to af and ac 
correspondingly. 

Let us consider now the "odderon"-Regge trajectory - the leading glu- 
onic trajectory with negative signature and C-parity. Mass of the lowest 3g 
glueball with spin 3 corresponding to this trajectory has been estimated 
above and found to be large 3.6 GeV in accord with lattice data. The 
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slope a' 3g for this trajectory should be equal to the one of gg-trajectory [] 
and thus the intercept of the nonperturbative glueball "odderon" is very low 
«3g(0) ~ —1.6. Mixing with gg-trajectories (u, (p) is much smaller than in the 
Pomeron there is no crossing of the odderon and (u, if) trajectories 

in the small t-region. 

The main results of the paper can be summarized as follows. We calcu- 
lated the 2g and 3g glueball spectrum analytically and compared resulting 
masses with lattice data, finding an excellent agreement, within 5% on av- 
erage for ratios of glueball masses. We stress that our spectrum contains 
no arbitrary fitting parameters, since all masses are expressed in terms of 
the string tension, which is fixed by the experimental meson Regge slope. A 
smallness of Coulomb-type contributions to glueball masses is pointed out. 
The spin splittings of glueball masses were obtained from first perturbative 
corrections calculated with nonperturbative wave functions. There is a good 
agreement on spin splittings between our calculations and lattice data as well. 
This indicates that the main ingredient of the glueball dynamics is the ad- 
joint string occurring between gluons in the two-gluon glueballs. The string 
dynamics implies that glueball spectrum corresponds to straight-line Regge 
trajectories, which have the Regge slope equal to | of that for meson tra- 
jectories. For the Pomeron Regge trajectory we found that with an account 
of only gluonic contributions (both nonperturbative and perturbative) the 
intercept is below unity and in order to obtain a phenomenologically except- 
able value of the intercept it is necessary to take into account mixing between 
gluons and qq pairs. We note that both nonperturbative (string dynamics, 
quark loops) and perturbative effects are important to obtain op(0) > 1. 
In this approach the "soft" and "hard" dynamics are strongly intermixed to 
produce the leading Pomeron pole. 

The partial support from the RFFI grant 96-15-96740 is gratefully ac- 
knowledged. One of the authors (A.K.) acknowledges a support of the NATO 
grant OUTR.LG 971390. 



3 The situation is analogous to the case of gg(meson) and qqq (baryon) Regge trajec- 
tories, since both for baryons and 3g glueballs at large L the structure q — 2q(g — 2g) is 
energetically preferrable. 
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Table 1 

Spin averaged glueball masses M G j \[d~j 



Quantum 
numbers 


This 

work 


Lattice data 


ref. [3] 


ref. [2] 


2 gluon 

states 


I = 0, n r = 


4.68 


4.66±0.14 


4.55±0.23 


I — 1, n r — 


6.0 


6.36±0.6 


6.1 ±0.5 


Z = 0, n r = 1 


7.0 


6.68±0.6 


6.45±0.5 


Z = 2, n r = 


7.0 


9.0 ±0.7(3++) 


7.7 ±0.4(3++) 


I — l,n r — 1 


8.0 




8.14 ±0.4(2*"+) 


3 gluon 

state 


K=0 


7.61 




8.19±0.48 



Table 2 

Comparison of predicted glueball masses with lattice data (for 

a f = 0.238GeV 2 





M(GeV) 
This work 


Lattice data 


M[G]/M[0"+] 


Difference 


ref. [2] 


ref. [3 


This work 


ref. [2] 


0++ 


1.58 


1.73±0.13 


1.74±0.05 


0.62 


0.67(2) 


-7% 


0++* 


2.71 


2.67±0.31 


3.14±0.10 


1.06 


1.03(7) 


3% 


2++ 


2.59 


2.40±0.15 


2.47±0.08 


1.01 


0.92(1) 


9% 


2++* 


3.73 


3.29±0.16 


3.21±0.35 








o-+ 


2.56 


2.59±0.17 


2.37±0.27 










3.77 


3.64±0.24 




1.47 


1.40(2) 


5% 


2-+ 


3.03 


3.1±0.18 


3.37±0.31 


1.18 


1.20(1) 


-1% 


2 -+* 


4.15 


3.89±0.23 




1.62 


1.50(2) 


8% 


3++ 


3.58 


3.69±0.22 


4.3±0.34 


1.40 


1.42(2) 


-2% 


1— 


3.49 


3.85±0.24 




1.36 


1.49(2) 


-8% 


2~ 


3.71 


3.93±0.23 




1.45 


1.52(2) 


-1% 


3~ 


4.03 


4.13±0.29 




1.57 


1.59(4) 
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Figure Caption 



Fig.l Vacuum trajectories before mixing (dotted lines) and after mixing 
(solid lines). 
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